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QCD is investigated at finite temperature using Wilson fermions in the fixed scale approach. A 
2+1 flavor stout and clover improved action is used at four lattice spacings allowing for control 
over discretization errors. The light quark masses in this first study are fixed to heavier than phys- 
ical values. The renormalized chiral condensate, quark number susceptibility and the Polyakov 
loop is measured and the results are compared with the staggered formulation in the fixed A^, 
approach. The Wilson results at the finest lattice spacing agree with the staggered results at the 
highest A^,. 
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1. Introduction 

Extensive experimental work is currently being done with heavy ion collisions to study the 
QCD transition, most recently at the Relativistic Heavy Ion Collider, RHIC and at the Large Hadron 
Collider, LHC. Both for the cosmological transition and for RHIC/LHC, the net baryon densities 
are quite small, thus the baryonic chemical potentials (/i) are much less than the typical hadron 
masses, jJ. is below 50 MeV at RHIC, even smaller at LHC and negligible in the early universe. 
Thus in this study we stay at /i = 0; for a review on thermodynamics see [1]. 

When one analyzes the absolute scale or any other question related to the T >0 QCD transition 
for the physically relevant case two ingredients are quite important. 

First of all, one should use physical quai^k masses. The nature of the transition strongly de- 
pends on the quark mass. Lattice studies and effective models showed that in the three flavor 
theory for small or large quark masses the transition is a first order phase transition, whereas for 
intermediate quark masses it is an analytic crossover. 

Secondly, the nature and other characteristics of the T > QCD transition is known to suffer 
from discretization errors [2, 3]. Let us mention one example which underUnes the importance 
of removing these discretization effects by performing a controlled continuum extrapolation. The 
three flavor theory with a large, a w 0.3 fm lattice spacing and standard staggered action pre- 
dicts a critical pseudoscalar mass of about 300 MeV [4]. This point separates the first order and 
cross-over regions. If we took another discretization, with another discretization error, the critical 
pseudoscalar mass turns out to be much smaller, well below the physical pion mass of 135 MeV. 
The only way to determine the physical features of the transition is to carry out a careful continuum 
limit analysis. It can be safely done only in the so-called scaling regime. To carry out a controlled 
continuum extrapolation at least three lattice spacings in the scaling regime are needed, because two 
points will always lie on a 2-parameter curve describing the conections to the continuum results. 

Numerically it is very demanding to fulfill both conditions. There are only a few cases for 
which this has been achieved. However it is important to note that fulfilling the second condition 
(at least 3 lattice spacings) without fulfilling the first one (physical quark masses) still leads to 
universal results. In other words continuum extrapolated results with non-physical quark masses 
are universal. These results are not the same as they are for physical quark masses, but they are well 
defined and unique. Contrary to this universality, fulfilling the first condition (physical quark mass) 
but not the second one (at least 3 lattice spacings) leads to non-universal, non-physical results. 
These results still have unknown discretization errors. 

For this reason in this first study of QCD thermodynamics with Wilson fermions we have 
chosen to work at heavier than physical quark masses but simulate at four lattice spacings in order 
to approach the continuum limit. 

In this paper we determine the temperature dependence of the chiral condensate, strange quark 
number susceptibility and Polyakov-loop in 2-1-1 flavor QCD. We use Wilson fermions with six 
stout smearing and tree-level clover improvement in the quark sector and tree-level improved fields 
in the gauge sector; for the details of the action see [5]. 

The structure of the paper is as follows. After this brief introduction the main features of the 
action and run parameters are listed in Section 2. Our choice of renormalization procedures for the 
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Table 1: Simulation parameters. The A', values used for the finite temperature runs and the values used for 
the zero temperature runs are separated by a comma. For finite temperature only even A', values were used. 
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Table 2: Mass ratios from spectroscopy. The lattice spacings are given by ma = 1672 MeV. 

various measured quantities are summarized in Section 3. The results are given in Section 4. In 
Section 5 we summarize and provide a short outlook. 

2. Simulation points and techniques 

The gauge action is the tree level improved Symanzik action [6], the fermionic action is a 2-1-1 
flavor stout [7] and clover [8] improved Wilson action. The number of smearing steps is six and 
the smearing parameter is set at 0. 1 1 ; for the first large scale simulation employing stout smearing 
see [11]. The clover improvement is tree level, i.e. = 1. In order to speed up the molecular dy- 
namics multiple time scales [9], Omelyan integrator [10] and even-odd preconditioning was used. 
For further details on the algorithm see [5] . The finite temperature simulations were performed in 
the fixed scale approach [12]. 

The simulation parameters are summarized in Table 1 . The simulations were performed at four 
lattice spacings and dimensionless ratios of hadron masses were used to define the line of constant 
physics. The mass of the Q. baryon, ma = \612MeV , was used to set the scale. The ratio of the 
pion and omega mass was kept at heavier than physical value, while the mostly light quark mass 
independent combination {2m\ — m\)/m^Q^ was kept close to its physical value. More precisely 
the combination which was tuned the most precisely is {2m\ — m1)/m1 which can also be used to 
define a physical ratio ms/mud while m^ itself was fixed by the requirement that as mud is lowered to 
the physical point the kaon mass becomes physical too. The large lattice volumes make sure finite 
volume effects are negligible, in all cases mj^L > 7.9. Table 2 summarizes the zero temperature 
hadron spectroscopy results. For more details see [13, 14, 15]. 

3. Renormalization 

Three quantities were measured, the chiral condensate {{yy)), the strange quark number sus- 
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ceptibility {Xs) and the Polyakov loop (P). Xs does not require renormalization but the others do. 

3.1 Chiral condensate 

In this section the additive and multiplicative renormalization of the chiral condensate is sum- 
marized briefly, for more details see [16] which is based on [17]. 

The additive divergence in the bare condensate is 0{a^^) which at finite temperature can be 
removed by subtracting the zero temperature value. The multiplicative divergence also needs to be 
removed by multiplying an appropriate expression of the quark masses. This step is complicated 
by the fact that Wilson fermions break chiral symmetry explicitly but nevertheless the axial Ward 
identity ensures that the following two combinations are finite and agree in the continuum limit 
[16], 

mR{f\i/)R{T) = mpcAcZA^w(.T) + 0{a) (3.1) 
mR{Y\i/)R{T) = 2NfmlcAcZl^PPiT) + 0{a) 

where Za is a finite renormalization constant, 

A^^(r) = {yoYo){T)-{Yo¥o){T = 0) (3.2) 
App{T) = I d^x{Po{x)Pa{0)){T)- j J^^(PoWPo(0))(r = 0) 

are the subtractions needed to cancel the additive divergences and {yqYo) and Pq{x) are the bare 
condensate and bare pseudo scalar density, respectively [16]. In this first study we did not measure 
Za but rather used the two expressions (3.1) to eliminate Za (and also mpcAc) and end up with the 
finite expression 

3.2 Polyakov loop 

The divergence of the bare free energy means that the Polyakov loop also needs to be renor- 
malized. A convenient choice of renormalization prescription is demanding a fixed value for the 
renormalized Polyakov loop at a fixed but arbitrary temperature > Tc. Then the renormalized 
Polyakov loop Pp is given by 

in terms of the bare Polyakov loop Po(P). We choose = 0.143ma and P* = 1.2. Other choices 
would simply correspond to other renormalization schemes [20] . 



4. Results 

We have measured three quantities at four lattice spacing. The renormalized chiral condensate 
is sensitive to the chiral transition whereas the renormalized Polyakov loop and the strange quark 
number susceptibility are sensitive to the remnant of the confinement-deconfinement transition. 
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Figure 1: Renormalized chiral condensate. The four lattice spacings with Wilson fermions (in the fixed 
scale approach) are compared with three lattice spacings with staggered fermions in the fixed A^, approach. 



Finite continuum results should agree between different discretizations provided the same 
renormalization conditions are used. We have imposed the same conditions on staggered simu- 
lations hence these can be compared with the Wilson results; for a recent review on staggered 
thermodynamics see [18]. 

Figure 1 shows the renormalized chiral condensate at the four different lattice spacings to- 
gether with staggered results at three fixed A^, values, 8, 10 and 12 (the staggered results were 
obtained using the fixed Nt approach, for more details see [19, 20, 21, 22]). As is apparent from 
the figure the two finer lattice spacing values are small enough to observe only a small change in 
the chiral condensate at any temperature. The rougest lattice spacing on the other hand is seen to 
include large discretization effects and is probably not in the scaling region. What is also clear is 
that results atNt <S are not trustworthy at any of the lattice spacings due to a too rough resolution 
of the imaginary time direction (temperature). At only moderately high and small temperatures, 
particularly around Tc, the three finest lattice spacings are probably in the scaling region and con- 
tinuum extrapolation is possible. 

The strange quark number susceptibility is a sum of two contributions, the connected and 
disconnected terms. The disconnected part is a very noisy quantity (as usual) and a large number 
of random vectors are needed in order to evaluate it precisely [23] . The results for the four lattice 
spacings are shown on Figure 2 together with the staggered results. We plot Xv/^^ normalized by 
its Stefan-Boltzman limit value. 

The third quantity we measured is the renormalized Polyakov loop which is sensitive to the 
confinement-deconfinement transition similarly to the quark number susceptibility. The results aie 
shown on Figure 3. 
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Figure 2: Strange quark number susceptibility. The four lattice spacings with Wilson fermions (in the fixed 
scale approach) are compared with three lattice spacings with staggered fermions in the fixed A^, approach. 
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Figure 3: Renormalized Polyakov loop. The four lattice spacings with Wilson fermions (in the fixed scale 
approach) are compared with three lattice spacings with staggered fermions in the fixed Nt approach. 
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5. Summary 

In this work we have started our study of 2+1 flavor QCD thermodynamics with Wilson 
fermions. Special emphasis is placed on obtaining results in the continuum limit, i.e. the sim- 
ulations are performed at four lattice spacings. In this first study the light quark masses were 
kept above the physical value but since continuum results are universal even at non-physical quark 
masses a comparison is possible with continuum staggered results. We found nice agreement be- 
tween the finest lattice spacings of the Wilson and staggered formulations. In the future we plan to 
extend our results to lower pion masses and also plan to use (3. 1) for (i/^i/'^) after measuring Za and 
fnpcAC because a priori it is not clear which formula has the better scaling property. 
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